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Abstract 

We present a practical and efficient means to compute the singular value decomposition (svd) of a 
quaternion matrix A based on bidiagonalization of A to a real bidiagonal matrix B using quaternionic 
Householder transformations. Computation of the SVD of B using an existing subroutine library such 
as LAPACK provides the singular values of A. The singular vectors of A are obtained trivially from the 
product of the Householder transformations and the real singular vectors of B. We show in the paper that 
left and right quaternionic Householder transformations are different because of the non-commutative 
multiplication of quaternions and we present formulae for computing the Householder vector and matrix 
in each case. 

1 Introduction 

The singular value decomposition (svd) of a quaternion matrix was first described theoretically in 1997 by 
Zhang 13, Theorem 7.2]. Zhang demonstrated the existence of the quaternion SVD using an isomorphism 
between a quaternion matrix and a complex matrix known as the complex adjoint. The details of the complex 
adjoint matrix need not concern us here, but it is important to note that this matrix is redundant — it 
contains twice the number of real values as the corresponding quaternion matrix. That is, every quaternion 
is represented by four complex numbers in the complex adjoint, rather than the two that are sufficient, 
the pattern of sign/conjugation and the position of these complex numbers preserving the properties of the 
quaternion, so that computing with the complex adjoint is equivalent to computing with the quaternion 
matrix. Therefore, computing the SVD of the complex adjoint matrix is equivalent to computing the SVD of 
the corresponding quaternion matrix, but since the SVD of the complex adjoint contains the singular values 
of the quaternion matrix twice, alternate singular values of the complex adjoint have to be discarded. Details 
on using the complex adjoint matrix to compute the quaternion SVD may be found in [H] but this method 
is effectively superseded by the results in this paper*. 

It has been noted by the authors of a quaternion QR algorithm 2 that using complex adjoint matrices 
to compute quaternion matrix decompositions can lead to loss of accuracy, since an algorithm applied to a 
complex adjoint does not necessarily preserve accurately the pattern of values in the complex matrix. For 
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this reason, and because the computation inherently involves twice as much work as a direct quaternion 
computation, we have sought direct quaternion methods to compute the quaternion SVD. 

In this paper, we present a method to compute the SVD of an arbitrary quaternion matrix A which is 
both elegant and fast, based on transformation of A to real bidiagonal form using quaternionic Householder 
transformations. We show in the paper that quaternionic Householder transformations have two forms, 
dependent on whether the Householder matrix is to be applied on the left or right of the vector to be 
transformed (this is because quaternion multiplication is not commutative). Using left and right Householder 
transformations, the bidiagonalization of a quaternionic matrix may be represented by: HAG = B, where 
H (G) is the product of the left (right) quaternion Householder matrices required to transform A to the 

real bidiagonal matrix B. Since H and G are unitary, we have also*: A = H BG . The SVD of the real 
bidiagonal matrix B may be computed using existing algorithms and code for real matrices. This is usually 
represented by B = USV T , where U and V are orthogonal real matrices. We show that the SVD of the 
quaternion matrix A is then: 

A = H T USV T G T (1) 

Thus the singular values of the quaternion matrix A are the same as the singular values of the bidiagonal 
real matrix B, and the singular vectors of A are obtained trivially from the products H U and V T G . 

The most significant aspect of this method is that the SVD of the real matrix can be computed using an 
existing subroutine or library. Considerable effort has been expended by others (for example, the authors 
of the lapack library £Q) to compute the SVD of real (and complex) matrices efficiently and accurately and 
the results of this effort can be exploited in computing quaternion SVDs without having to write quaternion 
versions of the SVD computation itself. The idea of computing the SVD using reduction to bidiagonal form is 
well known, but it is much less well known that the bidiagonal form may be real, even for a complex matrix*. 

The sequence of topics in the rest of the paper is as follows. In |2|we discuss the generalization of the 
Householder transformation to quaternion vectors; we prove that the transformation exists, and is unitary 
for both the left and right cases; and we present an algorithm for computing the quaternionic Householder 
vector in each case. In |3] we prove the existence of the real bidiagonalization of an arbitrary quaternion 
matrix using the quaternionic Householder transformations. In ^1 we give an algorithm for bidiagonalizing 
a quaternion matrix. We then show in §3 that the singular values of the bidiagonalized matrix are identical 
to the singular values of the original quaternion matrix, and we present an algorithm for computing the 
quaternion SVD using the quaternion bidiagonalization in 21 and any existing algorithm for the SVD of a real 
matrix. 

We mention here for completeness that it is also possible to compute the real bidiagonal matrix B using 
a quaternionic form of the Golub-Kahan-Lanczos bidiagonalization E] . This is conceptually simpler than 
the use of Householder transformations, since the only change needed to generalize the algorithm to handle 
quaternion matrices is the definition of the norm of a quaternion vector, which we give here in Lemma |21 
However, the Golub-Kahan-Lanczos algorithm is not attractive because of its computational complexity. 
For each non-zero element of the bidiagonal matrix B generated, the corresponding quaternion vector is 
generated, but to yield accurate results, these quaternion vectors must be re-orthogonalized against all the 
previously generated vectors, and this requires considerable computation. We have implemented the Golub- 
Kahan bidiagonalization for quaternion matrices and verified that it works, but the method presented in 
this paper (which is also a generalization of the Householder method first presented in j3j) is much faster, 
because it does not require re-orthogonalization. 

We assume familiarity with quaternions and quaternion matrices. The general concepts of quaternions 
are available in several books [71 1121 IH], The paper by Zhang JH] is currently the most convenient and 
comprehensive summary of material on quaternion matrices. Some elementary results which we need in the 
rest of the paper are given in Appendix 1X1 

tin this paper we use the overbar to denote a quaternion conjugate and a superscript T to denote a transpose. We avoid the 
standard superscript asterisk in our work, reserving it to denote a complex conjugate in work with complexified quaternions. 

t Interested readers may like to inspect the code of the LAPACK routine ZGEBRD or refer to the relevant section of the LAPACK 
Users Guide [Q available online at http://www.netlib.org/lapack/lug/node53.html 
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2 Quaternion Householder transformations 



The Householder transformation |S] is based on a unitary matrix of the form H = (i uu T ), where ||u|| = 
\f2. u is called a Householder vector. H is used to transform a vector, either by multiplication on the left by 
H for a column vector, or by multiplication on the right by H T for a row vector. Normally the vectors are 
parts of a row or column of a matrix that is transformed by a succession of Householder transformations. The 
modification usually consists of setting to zero multiple elements of the vector, while preserving the norm, 
but the general case is interesting. Given a unit vector v and a row or column vector a to be transformed, 
the transformation effectively generates a copy of v scaled by the norm of a, that is: Ha = ||a||v. Thus if v 
contains one non-zero element with a modulus of 1, the transformed vector will have one non-zero element in 
the corresponding position with a modulus equal to the norm of a and sign equal to the sign of the non-zero 
element of v. 

The Householder transformation was generalised to the complex case by Morrison JU| . He introduced a 
complex scalar with unit modulus, denoted £. The Householder matrix and vector then become complex: 
H = z (i — uu T ) where z = l/£. Now, if v has non-zero real elements, and a unit norm, then the transformed 
vector will be zero in all positions except those where v was non-zero, and the non-zero elements of the 
transformed vector will be real. This property follows from the fact that the transformation generates a 
copy of v scaled by the norm of a — if the elements of v are imaginary or complex, then so will be the 
corresponding elements of the transformed vector. Note that by using a real vector v the transformation 
results in a real modified vector. Therefore a sequence of Householder transformations may be used to 
bidiagonalize a complex matrix directly to a real bidiagonal matrix. 

In the complex case, z commutes with elements of the matrix I— uu T (and with elements of a). Therefore, 
the same matrix is used whether multiplied on the right or left (but as noted above, if multiplied on the 
right it must be transposed). 

In the quaternion case, we must distinguish between a left Householder transformation (to be applied 
to a column vector) and a right Householder transformation (to be applied to a row vector). To make this 
distinction clear, we use a different symbol in each case, H for a left Householder matrix, and G for a right 
Householder matrix. In computing the Householder matrix, we must be careful to compute and use z and u 
correctly, because they do not commute with elements of a. It is worth emphasising that the terms left and 
right are used in the matrix multiplication sense and not in the quaternion multiplication sense, although 
the difference between the two cases arises from the non-commutativity of quaternion multiplication. 

In this section, we generalize the Householder transformation to the case of an arbitrary quaternion 
vector a, but we are forced to restrict the vector v to be real in order to solve for the Householder vector 
u and the scalar £. This is not an important restriction, since in the paper we are specifically interested in 
the case where v and the transformed vector are real. We then show how the left and right Householder 
transformations are related. 

In the appendix, we present verification of several important properties of the Householder matrix in the 
quaternion case. Since these properties generalize from the complex case, their proofs are not vital here, but 
they are included in the paper for completeness. 

The authors of [2] used quaternion Householder transformations in their QR algorithm and gave an 
algorithm for computing a Householder vector, based on that of Morrison. However, they did not discuss 
the difference between left and right Householder transformations. 

2.1 Left Householder vector and matrix 

In this case the Householder matrix H is applied to the column vector a on the left, and the scalar z is 
multiplied on the left of the matrix I — uu . Thus H = z (I uu . 

Theorem 1. Given an arbitrary quaternion vector a £ H r , and a real vector v € W. r with unit norm; 
there exist a quaternion vector u G H r with norm ||u|| = v2 and a unit quaternion scalar z, such that 
Ha = z (i — uu T ) a = ||a||v. 
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Proof. The proof is by construction, and is based on that of Morrison ^U] for the complex case, with 
the necessary adjustments to allow for non-commutative multiplication. We use Morrison's symbols, and 
therefore we write £ = and a = ||a||. Therefore, we have: 

(I - uu T ) a = (va (2) 
Multiplying out the parentheses on the left, we obtain: 

a — uu T a = (va (3) 
and then multiplying from the left by a T , we obtain: 

rj-\ rj-i rrt rrt , . 

a a — a uu a = a Qva (4) 

Now, a T a = a 2 by Lemma El and a T uu T a is the product of a quaternion a T u with its conjugate. Denoting 
the quaternion^ u T a by fi, we then have: 

a 2 - \[i\ 2 = a T Cva (5) 

The left-hand side of this equation is real, and therefore so must be the right-hand side. Therefore the 
product a T £ is real, since v and a are both real. At this point, Morrison was able to reorder the terms and 
solve for £ by writing a T v as the product of a positive real (modulus) with a unit complex number. We 
cannot do this because C does not commute with the other terms. However, we can re-order terms by the 
quaternion conjugate rule (qTq2 — 92 <7i)- Thus: a T £v = £a T v. Since the result here is real, we may drop 
the overall conjugate and obtain^: a T £v = £a T v. Following Morrison, we now represent a T v as the product 
of a real modulus with a unit quaternion, and write it as a T v = rw, where r is the modulus, and uj is a unit 
quaternion. Then, the fact that £a T v is real implies that ( = ±w since £ has unit modulus. This implies 
that £ = ±w. Substituting this result into equation \5\ and remembering that lj is a unit quaternion, we 
obtain: 

a 2 — |/i| 2 = ±Uruja — ± ar (6) 

which gives on re-arrangement: 

H 2 = « 2 T ar (7) 

We now make an arbitrary choice, and choose fi to be real (we could have done this earlier), hence we have: 
/i = -J a(a =F r). Of the two possibilities here, we follow Morrison and choose to add the terms inside the 
parentheses (both are positive, and if we subtract them we might get a small result in some cases). Therefore 
we have fi = W a(a + r) and the choice we have just made also implies that Q = —lo. Finally from equation 
13 we can obtain u: 

a — (va — uu T a = u/i (8) 

and remembering that we have chosen /i to be real: 

u = - (a - Cva) (9) 
A* 

The sequence of steps necessary to compute u and £ is therefore as follows: 

a=|NI (10) 



§In fact, we make the choice later that fi is real. 

^It is at this point that we are forced to impose the restriction that v be real, since otherwise we could not re-order the 
right-hand side to put £ on the left of a T . 
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If a = choose u to be a zero vector and £ = 1 (this causes the Householder matrix H to be an identity 
matrix), otherwise continue: 



r = 


|a T v| 






c = 


I a v 




I r 


M = 


y/ a(a + r) 






u = 


— (a — (va) 



r = 
r > 



(11) 
(12) 

(13) 
(14) 

M 

It remains to show that ||u|| = y2- From LemmaHJ ||u|| 2 = u T u, hence: 

||u|| 2 = -(a - Cva) T - (a - Cva) = \ (a T - Cv T a) (a - (va) (15) 
MM M 

Multiplying out the terms in parentheses, and re-ordering terms that commute: 

||u|| 2 = -j (a T a-Cv T aa-a T vCa + CCv T va 2 ) (16) 
M 

a T a = a 2 by Lemma|21 CC — IC| 2 — 1 an d vTv = ll v l| 2 — 1; an d M 2 = a ( a + r )j hence: 

2 a 2 — £v T aa — a T v(a + a 2 2a — £v T a — a T v£ 
a(a + r) a + r 

Since v is real, v T a — a T v, and we have: 

u 2 = ^— i 18 

a + r 

Finally, from eauation ll2l a T v = — r(, and therefore a T v = — r£ since v and r are real. Therefore: 

2 a + rCC + rCC = 2a ± 2r =2 
a + r a + r 

as required. □ 
2.2 Right Householder vector and matrix 

In this case the Householder matrix G is applied to the row vector a T on the right, and the scalar z is 
multiplied on the right of the matrix I — uu T . Thus G = (I — uu T ) z. 

Theorem 2. Given an arbitrary quaternion vector a T S IF, and a real vector v T £ R r with unit norm, 
there exist a quaternion vector u T € HT with norm ||u|| = \/2 and a unit quaternion scalar z, such that 
a T G = a T (I - uu T ) z = ||a||v T . 

Proof. This theorem follows from Theorem ^ 

a T (l-uu T )z = |M|v T 

Place z on the left using the quaternion conjugate rule (<?i<72 = §2 Qi) '■ 



za T (I - uu ) = llallv 



T 
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Transpose both sides: 



za T (I — uu T ) = ||a||v 



Using Lemma m we place a on the right: 



z (I — uu T ) a = ||a||v 
Since z is a scalar we may take it under the overall conjugate transpose: 

T 

z (I — uu T ) a=||a||v 

We recognise the left-hand side to be H a as in Theorem ^ That is, we may compute a right Householder 
matrix by taking the conjugate transpose of the left Householder matrix computed from a conjugated (and 
transposed) vector a. □ 

For the sake of completeness, we give without proof the algorithm required to compute a right Householder 
transformation directly. The sequence of steps necessary to compute u T and £ is as follows: 

a=||a|| (20) 

If a = choose u to be a zero vector and £ = 1, otherwise continue: 

r = |v T a| (21) 

r ' r = ° (22) 
l_v_a . r>Q \ J 



/j = yja(a + r) (23) 

u T = i(av T C-a T ) (24) 
H 

3 Existence of the bidiagonalization of a quaternion matrix 

Golub and Kahan in their 1965 paper |3j demonstrated the existence of the bidiagonalization of a real 
matrix and showed how to decompose an arbitrary matrix into the product of two orthogonal matrices and a 
bidiagonal matrix, either using the method now known as the Golub-Kahan method, or by using a sequence 
of Householder transformations. We show here, using quaternion Householder transformations, that the same 
decomposition is valid for an arbitrary quaternion matrix and, more significantly, that the bidiagonal matrix 
may be real. This latter property follows from the behaviour of the quaternion Householder transformations 
in fJ21when the vector v is real. 

Theorem 3. Given an arbitrary quaternion matrix A £ H rxc with r rows and c columns, there exists a pair 
of unitary quaternion matrices L £ M rxr and R £ H cxc , and a real bidiagonal matrix B £ R rxc such that 
LAR = B. B is upper or lower bidiagonal dependent on the relative magnitudes of r and c. Specifically, if 
c> r (r > c), B is non-zero only on its diagonal and super-diagonal (sub-diagonal respectively). If r = c, B 
may be non-zero on either its sub- or super- diagonal arbitrarily. 

Proof. The proof is by induction and construction. We show that the first column of A may be transformed 
so that the first element is the only non-zero element and is real. We then show that the remainder of the 
first row may be similarly transformed by transposing and conjugating the remaining columns and applying 
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Figure 1: Construction of an upper bidiagonal matrix from columns and rows each with one non-zero leading 
element. (Non-zero elements are marked by x.) 



the same process to the first column of the result. The process continues until the matrix to be transformed 
has one column only. Transforming this column completes the process of bidiagonalization. 

Assume for the moment that c > r and therefore that B is to be upper bidiagonal. Figure ^ shows 
how the upper bidiagonal matrix B may be constructed recursively from columns and rows with only one 
non-zero element. We may transform A using a left Householder transformation as defined in Theorem^^ 
that all elements of the first column become zero, except for the first element which becomes real (that is, it 
has zero vector part). The required left Householder matrix H is computed by choosing a equal to the first 
column of A, and v = (1, 0, 0, 0, ... , 0) T of length r, the same length as a. The left unitary matrix L is then 
equal to H. Multiplying A on the left by L transforms the first column (and modifies all the subsequent 
columns). Since v is real and has only one non-zero element, the transformed first column of A will also be 
real with one non-zero element as shown in Figure ^ The right unitary matrix R is an identity matrix. 

If A has only one column, it is now in (degenerate) real bidiagonal form as is evident from the first 
column of Figure ^ 

If A has more than one column, let A' be the matrix with r rows and c — 1 columns obtained by deleting 

/■ 

the first column of A. A' may be bidiagonalized by applying the above process to A' by LemmaQ] The 

T T — T T 

result will be: R' A' L' = B' where L' will be the same size as L, and R' will be one row and column 
smaller than R. Replacing A' by B' and multiplying L on the left by L' and replacing the submatrix 
R(2 — > c, 2 — > c) by R' completes the proof. 

In the case when r > c, we have from Lemma ^ 

— T — T — T — T 

R A L =B 

and thus L, R and B exist because we may form A . Since B is real, B = B T 

In the case where A is square, and r = c, we may choose B to be upper or lower bidiagonal arbitrarily. 
The simplest choice is upper bidiagonal. □ 

A dual proof is possible using right Householder transformations applied to rows of A. In this case the 

result would be lower bidiagonal, and the upper bidiagonal case would be handled by applying the process 

-r-T 
to A . 

We note, incidentally, that since real and complex numbers may be treated as subsets of the quaternions, 
the above proof is also valid for real and complex matrices (of course, the result is well-known in these cases). 
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4 Bidiagonalization using Householder transformations 



We give here a concise recursive algorithm to bidiagonalize an arbitrary quaternion matrix based on left 
transformations only in which rows are transformed by conjugating and transposing the matrix. Algorithm 
^is a reference algorithm based on explicit Householder matrices. It is easily implemented (we have done so 
in Matlab, as part of our quaternion library QD), but for large matrices, or matrices in which one dimension 
is large, the explicit Householder matrix may be too large to be computed explicitly. 

5 Quaternion Singular Value Decomposition 

The existence of the bidiagonalization of a quaternion matrix is not in itself sufficient to prove that computing 
the singular values of the bidiagonal matrix yields the singular values of the original quaternion matrix. We 
therefore prove that the singular values of the bidiagonal real matrix are the same as the singular values of 
the original quaternion matrix. 

Theorem 4. Given an arbitrary quaternion matrix A € H rxc ; and a real bidiagonal matrix B € M. rxc as 
defined in Theorem^ the singular values of A. are the same as the singular values o/B. 

Proof. From Theorem [21 there exist unitary quaternion matrices L and R that will transform A to B, that is 
LAR = B and since L and R are unitary, A = L BR . The singular value decomposition of B = USV T 

rp rp 

where U and V are orthogonal, hence A = L USV T R . From the uniqueness of the singular values, and 

rp rp 

from the fact that L U is unitary and V T R is unitary, it follows that £ contains the singular values of 
the quaternion matrix A. □ 

Algorithm |21 presents formally the steps needed to compute the quaternion singular value decomposition 
using the bidiagonalization presented in 21 The SVD of the real matrix resulting from the algorithm presented 
in 21 may be computed using any existing real SVD algorithm or routine, for example, the SVD function in 
Matlab, or a suitable lapack subroutine (note that there are specialised routines in LAPACK for computing 
the SVD of real bidiagonal matrices, and one of these should be used in preference to an SVD routine for 
general real matrices to exploit the structure of the bidiagonal matrix). 

6 Conclusions 

We have presented a practical and efficient method for computing the singular value decomposition of an 
arbitrary quaternion matrix A. This method bidiagonalizes A to a real bidiagonal matrix B using quaternion 
Householder transformations and a simple recursive algorithm based on that of Golub and Kahan [3]. The 
SVD of B, computed with any available SVD algorithm for real matrices, followed by the multiplication of 
each of two quaternion matrices by a real matrix, yields the SVD of A. 

The few previous methods for computing the quaternion SVD were based on complex adjoint matrices, 
or direct computation of the quaternion SVD itself. The former suffers from the serious disadvantage that 
the computation is not accurate, and the singular values and vectors are computed twice over. The latter 
suffers from the disadvantage that it is complex to code, and requires a deep understanding of the SVD. To 
achieve the same computational efficiency and accuracy as the best available real and complex SVD routines 
would require considerable effort in coding and testing which the method presented in this paper neatly 
circumvents. In its place, it is only necessary to code an efficient version of algorithms ^ an d 121 

A Appendix 

We gather here some elementary results needed in the paper. Most of these results would be trivial in the 
complex case, but because of non-commutative quaternion multiplication we must be careful. 
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Algorithm 1: Quaternion bidiagonalization 



A G H rxc 

L G H rXr , Bel rxc , r t gH cxc 



Input 
Output 

if c < r 

Bidiagonalize (A , L, B, R) 
else 

— T 

Bidiagonalize (A , R, B, L); 
end 

procedure Bidiagonalize ( Input 

Output 



B = B T 

: Aei™ 

: Lei rxr , BeM rxc , ReH cxc ) 



Transform the first column of A, where v is the vector (1, 0, 0, 0, . . .) T 

L = fT(A(l->r,l),v), A = LA 
R = I cxc 

If there is more than one column, transform the remainder of the 
first row by applying the procedure recursively to the conjugate 
transpose of the remaining columns: 

if Ol 

Bidiagonalize (A(l r, 2 -> cf , R' , T, I/) 
L = L'L 

A(l -> r,2 c) = T T 
R(2 -> c, 2 -> c) = R' 
end if 

B = Real part of A (The vector part of A will be zero.) 

function H(a£i m , vel m ) 

return H G H roXm as defined in using a and v. 

end 



Algorithm 2: Quaternion singular value decomposition 

Input : AGi™ 

Output : Uei rxr , SeM rxc , VeH cxc 



Bidiagonalize A using Algorithm Q to obtain L, B and R , such that 
LBR T = A 

Compute the SVD of B , to obtain Wei™, S, XeM cxc , such that 
B = WSX T 

U = L W V =X T R 
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Elementary operations on quaternion matrices have been well-studied and gathered together in Zhang's 
1997 paper |13| . These include the concepts of unitary and Hermitian matrices which generalize from 
complex matrices to quaternion matrices without difficulty. The product of unitary quaternion matrices is 
itself unitary. 

T — T — T 

Lemma 1. For any two quaternion matrices of compatible dimensions: AB = B A . 

Proof. This result was given by Zhang |131 Theorem 4.1 (2)]. □ 

Lemma 2. Given an arbitrary quaternion vector u, its norm is given by vtFu or by \J u T u. 

Proof. The norm of a vector ||u|| = \ u i\ 2 - The inner product of u with its transpose conjugate yields 

the sum of the elements of u with their conjugates. Since a quaternion multiplied by its conjugate gives the 
square of its modulus, the inner product gives the sum of the moduli squared, and taking the square root 
gives the norm of the vector. □ 

Lemma 3. Given an arbitrary quaternion vector u, the matrix uu T is Hermitian. 

Proof. Elements of the matrix on the diagonal consist of quaternions multiplied by their conjugates UiUi, 
and are therefore real. Elements not on the diagonal are quaternion conjugate pairs, since they are of the 
form either uiuf or uluj. □ 

Theorem 5. A quaternion Householder matrix, H, of the form I — uu T , is unitary if u is a quaternion 
vector with j|u|| = \/2~. 

Proof. It is sufficient to show that H 2 = I, since H = H from Lemma and subtracting an Hermitian 
matrix from an identity matrix gives an Hermitian matrix (self-evident): 



Theorem 6. A unitary quaternion matrix U scaled by a quaternion z with unit modulus is unitary whether 
z is multiplied on the left or right o/U. 

Proof. Consider first the case where z is multiplied on the left of U: 



H 2 = (I uu T ) 2 = I 2 uu T I Iuu T + uu T uu T = I - 2uu T + u (u T u) u T (25) 



The term in the parentheses on the right is the square of the norm of u by Lemma |21 and since ||u|| = \/2, 
we have: 




□ 




■T 



(27) 



and applying the quaternion conjugate rule (<7i<Z2 = <ii Qi)'- 




Clearly, then: 



(zljf (zU) = V T zzU = U T |z| 2 U = U T U = I 



(28) 
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and 



(zU) (zU) T = zUU T z = zlz = |z| 2 I = I (29) 



If z is on the right of the matrix U, modification of equation 1271 shows easily that (Uz) = zU and the 
results of multiplying the scaled matrix by its transpose conjugate are as in equations 1281 and 1291 □ 
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